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Abstract 



Adopting a procedure previously proposed to quantitatively study two- 
dimensional pion interferometry pj, an equivalent 2-D \ 2 analysis was per- 
formed to test the resolving power of that method when applied to less favor- 
able conditions, i.e., if no significant contribution from long lived resonances 
is expected, as in kaon interferometry. For that purpose, use is made of the 
preliminary E859 K + K + interferometry data from Si + Au collisions at 
14.6 AGeV/c. As expected, less sensitivity is achieved in the present case, 
although it still is possible to distinguish two distinct decoupling geometries. 
The present analysis seems to favor scenarios with no resonance formation at 
the AGS energy range, if the preliminary K + K + data are confirmed. The 
possible compatibility of data with zero decoupling proper time interval, con- 
jectured by the 3-D experimental analysis 0], is also investigated and is ruled 
out when considering more realistic dynamical models with expanding sources. 
These results, however, clearly evidence the important influence of the time 
emission interval on the source effective transverse dimensions. Furthermore, 
they strongly emphasize that the static Gaussian parameterization, commonly 
used to fit data, cannot be trusted under more realistic conditions, leading to 
distorted or even wrong interpretation of the source parameters! 
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I. INTRODUCTION 



The second-order interferometry of identical particles is a powerful tool for probing the 
space-time zone from which they were emitted 0. Almost two decades ago, it was suggested 
as a possible signature of the Quark Gluon Plasma (QGP), expected to be formed in high 
energy Nucleus-Nucleus collisions, by probing the expected large space-time dimensions of 
the emitting system at freeze-out. About ten years ago, when the first O + Au runs from 
CERN/SPS became available, there were expectations that we could be seeing its formation, 
particularly from tttt interferometry [4j]. However, due to limited statistics, the correlation 
function at that time had to be projected in one dimension only, leading to ambiguity 
in describing the overall behavior of interferometry data, i.e., they could be equally well 
described by two very distinct freeze-out scenarios ||. One of them reflected a dynamical 
model in which the pions were formed after the hadronization of the QGP and the other 
one considered, instead,the formation of a hadronic gas of resonances. 

On the other hand, several studies ||-§| have shown that dynamical models considering 
expanding systems, would lead to effects that could dramatically distort the two-particle 
correlation function. Among them, the most significant effect |J was caused by long lived 
resonances, which later decayed into the observed particles. As a side consequence of this 
study, it was suggested to use pion interferometry to probe resonance formation at ener- 
gies where their fractions were unknown ||. Once again ambiguity in separating different 
scenarios emerged, evidencing symptoms of urgency for very accurate and high statistics 
data, which has become available more recently, allowing for multi-dimensional analyses. 
Nevertheless, together with improved data, more precise theoretical and phenomenological 
tests were required, leading to the method suggested in Ref . |l] , in which a two-dimensional 
X 2 analysis was proposed to study the resolving power of pion interferometry. For that pur- 
pose, two dynamical scenarios were considered which predicted similar correlation functions, 
although the underlying decoupling geometries differed considerably. In one, long lived res- 
onances were neglected, while in the other, a resonance gas with fractions predicted by the 



Lund model |I| was considered. 

Nonetheless, to quantify the differences in terms of a \ 2 interferometric analysis, the 
contribution of long lived resonances decaying into pions seemed to be essential. This fact 
led to the question whether the resolving power of the method would remain high under 
less favorable conditions, i.e., if only shorter lived resonances would contribute to the par- 
ticle yield, as is the case of K + K + interferometry, This is precisely the goal of this paper. 
Furthermore, the method is applied to test the hypothesis of zero time emission interval, 
suggested by the experimental fit using 3-D static Gaussian parameterization 0. In this 
study the influence of the time emission interval on the transverse radius parameter emerged 
naturally and another very important point was clearly emphasized, i.e., the static Gaus- 
sian parameterization, popularly used to fit data, is usually misleading in more realistic 
situations, and results in distorted or even wrong interpretation of the source parameters! 
Prior to reach these points, however, we present a brief summary of the theoretical model 
underlying the analysis, the so-called Covariant Current Ensemble Formalism ||[1TJ and a 
brief review of the method discussed in Ref. This is then adapted to the present case, 
in which use is made of the preliminary E859 bidimensional data on K + K + interferometry 
from the AGS/BNL. 
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II. THE CO VARIANT CURRENT ENSEMBLE FORMALISM 



Under idealized conditions the correlation function, (^(fci, A^), of two identical bosons 
probes their decoupling or freeze-out space-time distribution, p(x), through C 2 (fci,fc 2 ) = 
l+\p(ki— k 2 )\ 2 - However, in actual high energy reactions, final state interactions, correlations 
between coordinate and momentum variables, and resonance production distort this ideal 
interference pattern, corresponding only to Bose-Einstein symmetry(see e.g. ||- |§). This 
may lead to erroneous interpretation about the underlying information on the decoupling 
geometry coming from the second-order interference pattern. In realistic cases, then, it 
is mandatory to employ more general formalisms flexible enough to include such 



non-ideal effects, reflecting model dependent scenarios. In the Covariant Current Ensemble 
formalism, the correlation function can be expressed as p.|TT 



where T(q) = (q c /q)/(e Qc ^ q — 1) is the Gamow factor that distorts the interference pattern 
due to final state Coulomb interactions, with q c = 2-Kam and q = {—(k\ — k 2 ) 2 ) 1 ^ 2 . 

In general, when resonances are produced, the complex amplitude, G(k 1 ,k 2 ), can be 
written as 

G(k u k 2 ) « (£ f(K + /r) (1 - tqUr/Tr)- 1 j^k^u^)) , (2.2) 

r 

where f(K + /r) is the fraction of the observed A" + 's arising from the decay of a resonance of 
type r, which freezes-out with final four velocity v%. It should be noted that, in the absence 
of resonances, the sum in eq.( |2.2|) reduces to only one term, f(K + ) = 1. The currents, 
jo{ufki), contain information about the production dynamics. 

The ensemble average in the above notation is performed by using the following param- 
eterization for the implicit break-up distribution |||6| 

D(x,p) <x exp{-£L _ (igpl _ Clztf - j} S (E - E P )6>M , (2.3) 

where r = (t 2 — z 2 )^ is the freeze-out proper time, and 77 = |log((t + z)/(t — z)),y = 
\ \og((E + p z ) f(E — p z )) are the space-time and momentum rapidity variables, respectively. 
The correlation between these rapidities is estimated from the Lund model to be Arj ~ 0.8 
0, Y c = 0.7,and y* m = 0. As regarding resonance fractions, the Lund model [[LCj in the 
AGS range suggests essentially two contributions for that scenario, i.e., that f{K^ irec1 ) = 0.5, 
f(K + /K*) = 0.5. 

We recall that the transverse momentum in the more general model proposed in Ref. 
is assumed to arise entirely from the finite momentum spread Ap of the pion wave-packets. 
It should also be clarified that this model coincides with the Covariant Current Ensemble 
formalism in the case of minimum packets, when associating the momentum spread to the 
so-called pseudo-temperature, T PS , through Ap 2 /m = T PS . This pseudo-thermal ansatz, 
however, was previously used in order to derive an analytical form for the correlation function 



11]. In the present analysis we are basically considering the Covariant Current Ensemble 
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formalism but, since numerical calculations are carried out from the start, we consider the 
full thermal ansatz instead, in which T is the effective inverse transverse mass slope from the 
experimental fit, i.e., T = 0.18 GeV |12|], corresponding to an average momentum (fcy) 



0.49 GeV/c. It should be added that no clear difference could be seen when comparing 
the correlation functions corresponding to the thermal versus pseudo-thermal cases, in the 
same kinematical region. The currents in the thermal model may be written covariantly as 



jo(k) = yjv^k^e' 



By earring out the ensemble average in Eq. ( |2.2| ) with the aid of ( |2.3| ) and of jo(k) 



defined above, we obtain the expressions for G(ki, k 2 ) used in the numerical calculations: 

2 d2 /a r°° t 2 r+°° _(y-v*) 2 r+oo (v-y) 2 . . ... 

G(k 1 , k 2 ) oc e - q T R T/ A / Tdre-— 2 \ dye \ dr) e ~^r e M<iocosh v - qLSi nh v ) 

JO J —00 J— 00 

f{K + /r) (1 - iqu r /T T y l yj ' [m 1T cosh(y r - y 1 )}[m 2T cosh(y r - y 2 ] 

r 

exp{—miTCOsh(y r — yi)/2T} exp{—m2TCOsh(y r — y 2 )/2T} . (2.4) 

The single inclusive kaon distribution in this notation is P\{ki) oc G{k il hi), which can be 
written, with the help of eq. (|2.3|) , as 

roo T 2 r+oo (y-y*) 2 r+oo (v-y) 2 

G(ki,ki)(x / t dre'^ 1 / dye 2Y ? / dr] e 2A ^ 2 

JO J— 00 J— 00 

Yf(K + /r)[m iT cosh(y r - yi)} exp{-m iT cosh(y r -y^/T} . (2.5) 

r 

The aim of the present study is to test if multidimensional kaon interferometry can 
discriminate scenarios including resonances from those in which they are absent, even in the 
much less striking limit of no significant long lived resonance contribution to the kaon yield. 
For doing this, we apply the method suggested in Ref. |TJ to extract the rms transverse 
radius, Rt, at decoupling and the rms decoupling proper time interval, At. Note that we 
assume implicitly that the chaoticity parameter A = 1 throughout our analysis. 



III. x 2 ANALYSIS 

To compare theoretical correlation functions with data projected onto two of the six 
dimensions, we must compute the projected correlation function as 

r f I d 3 k 1 d 3 k 2 P 2 (k 1 , fc 2 ) A 2 (q T , q L ; k u k 2 ) 

<W«r,?Lj - j d s kld s k2 p l{kl) p l{k2) A 2 (q T ,q L ; h,k 2 ) ' 

where Pi and P 2 are the one and two kaon inclusive distributions, and A 2 is the experimental 
two kaon binning and acceptance function. All calculation were performed using the Monte 
Carlo importance sampling method adopted in the CERES code ||. 

The acceptance function for the E859 experiment was approximated || by 

A 2 (q T ,q L ;k u k 2 ) = A 1 (k 1 )A 1 (k 2 )Q(22 - \<pi - <p 2 \)5(q L - \k zl - k z2 \)5(q T - \k Tl - k T A) . 

(3.2) 



4 



FIGURES 






l(t 


60000 




50000 




40000 


; 9- 




- \ 




; z 


30000 


: ~U 


20000 




10000 





-40 



-20 



20 



40 



FIG. 1. Phase space generated by the CERES code, with the simplified cuts shown 
in eq. and eq.flOl). 



The angles are measured in degrees and the momenta in GeV/c. The single inclusive 
distribution cuts are specified by 

Ax(k) = 0(14 < 9 lab < 28)Q(p lab < 2.9 GeV/c)e(y min > 0.75) . (3.3) 

It can be seen from Figure 1 that the phase-space generated with the above cuts re- 
produces very closely that covered by the experiment [[J. Only a few excess of generated 
particles can be seen at low transverse momentum px (or fcjv). 

To assess the statistical significance of the differences between the fits obtained assuming 
resonance and non-resonance dynamics, we computed the x 2 goodness of fit, estimating this 
variable as previously 

X 1 ' J) {[AA(i, j)] 2 + [A/;- 1 ^, j)AB(i, j)] 2 } ' 1 ' 1 

where M x is a normalization factor which is chosen to minimize the average \ 2 an d de- 
pends on the range in the (qT,qi) plane under analysis. The indices i,j refer to the the 
corresponding qj-, qi bins, in each of which the experimental correlation function is given by 



(&A(i, 3 )\ 2 , fAB(i,j) 



(3.5) 
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The numerator A(i,j) ± AA(i,j) and denominator B(i,j) ± AB(i,j) in Eq. (|3.4j ) and 
( p.5[) were obtained from V. Cianciolo @>@) understanding that the data in this form 
are preliminary and subject to further final analysis. Use is made of its preliminary form 
mainly for testing the sharpness of the method. Note that in the present analysis we are 
not including the errors associated to the theoretical correlation function generated by the 
Monte Carlo importance sampling in CERES. All calculations, however, were performed by 
taking high statistics runs only, making it reasonable to consider those errors as negligible. 

Analogously to the procedure adopted in Ref. minimization of the average y 2 was 
performed by exploring the parameter space of Rt and At and computing the (x 2 )? averaging 
over a grid of nearly 30x30 bins in the (q?, qi) plane in the relative momentum region 
0.005 < qr,qi < 0.605 GeV/c, binned with 5qx = Sqi = 0.02 GeV/c. A very meticulous 
investigation was performed to find the most probable region where the minimum (Rt , Atq) 
could be located. 




FIG. 2. Zone in the (Rt,At) plane investigated, leading to the determination of 
the most probable region where the minimum (x 2 )? associated to (Rt,At), could be 
located. Part (b) corresponds to the dynamics ignoring the contribution of K* and 
part (d) to including their contribution to the kaon yield. Similarly, parts (a) and (c) 
correspond to non-resonance and resonance cases, respectively, but fixing At = 0, and 
optimizing only Rt. 
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Figure 2 illustrates most of the investigated region in the (Rt, At) plane. Parts (a) and 
(c) will be discussed latter. In the vicinity of the minimum, the parameters of the quadratic 
surface {x 2 (Rt, At)) = xLin + a (R-T ~ Rtq) 2 + /3(At — Ar ) 2 were determined. The results 
of such investigation are given in Table 1, where the radius parameters are measured in fm 
and time intervals in fm/c. 



TABLE 1: 2D-^ 2 Analysis of Kaon Decoupling Geometry 





No Res. (f Kdir = 1) 


LUND Res. {f Kdir = fx/K* = 0.5) 


Optimized Rt and At 


(Xmm)30x30 


1.03 


1.02 


(Xmin)lOxlO 


1.17 


1.30 


Rto 


2.19± 0.76 


1.95 ± 0.89 


Ar 


4.4± 2.0 


4.4± 2.6 


a 


0.0410 


0.0299 


P 


0.0058 


0.0034 


Optimized R T (At = 0) 


(Xmin) 30x30 


1.29 


1.33 


(Xmm)lOxlO 


4.04 


2.92 


Rq t 


~ 10.6 


~ 4.8 


a 


0.0003 


0.0280 



The errors appearing in Table 1 were estimated following the prescription of Ref. [|TJ, 
which considered the x 2 over N bins as a random variable and, for large N, approximated 
the distribution of the mean x 2 P er bin by P(x 2 ) exp[— (x 2 — l) 2 /2a 2 }, with rms width 



a = y2/N m 0.048, for the iV = 855 grid under consideration (i.e., subtracting from the 
original 900 the empty bins and the number of degrees of freedom consumed in the x 2 
analysis itself). Inserting the expression for (x 2 {Rt, At)) in the above paraboloid into the 
asymptotic form of the x 2 distribution for large N, the likelihood for the parameter Rt to 
have a value near the minimum is approximately oc exp[— o 2 (Rt — Rt ) a /2<J 2 ]. Therefore 
the estimated error on the radius is AR w { v / 2[r(3/4)/r(l/4)]a/a} 1 / 2 w 0.7( ( r/a) 1 / 2 , and 
similarly the error on the proper time interval is 0.7(er/ P) 1 ^ 2 . 

Comparing Table 1 with Ref. Jl|, we may see that the optimized value of At, the de- 
coupling time interval, is estimated to be about the same as in the pion case. However, 
the transverse size of the kaon emission region seems to be half that of the pions. This 
results agrees with the experimental fit to the data and, as was stated in Ref. 0, it could 
be reinforcing the suggestion in Ref. IfLil , according to which kaons could decouple earlier 
than pions due their small cross section for interacting with nuclear matter. 

We see from Table 1 that the optimization in both scenarios result in similar values for 
(x 2 ) over 855 bins, although using the optimized parameters, we see smaller (x 2 ) for the 
non-resonance scenario in a smaller (10 x 10) grid. Just to illustrate the similarities, we can 
see in Figure 3 the two-dimension correlation functions C(qT, <?l), corresponding to data and 
to the theoretical values generated with the optimized values shown in Table 1. 
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FIG. 3. The preliminary E859 Gamow corrected data are shown in part (a). Part 
(b) shows theoretical correlation functions, C{qr, filtered with the E859 acceptance 
for the case with no resonances (NR) and part (c) corresponds to the inclusion of Lund 
resonance (LR) fractions; the corresponding distribution of X 2 (qt,Ql) are m (f) and 
(g), respectively. Similarly, when fixing At = 0, results for the generated C{qT,qL) 
are shown in parts (d) (NR) and (e) (LR), with X 2 (QT,qL) distribution in (h) and (i), 
respectively. 

From the above discussion, similarly to what happened in the pion case, we see that not 
enough separation is found, neither from the 2-D projection alone, nor by conjugating it 
to the average \ 2 analysis. However, in Ref. JIJ it has already been recalled that a most 
direct measure of the goodness of fit could be achieved by means of n a = \(Xmin) ~ the 
number of standard deviations from unity of the average x 2 P er degree of freedom. Since no- 
depends on the range of q under analysis, we followed the steps of Ref. |TJ and studied its 
behavior by varying the range of the analysis to restricted (q?, ?l) domain, ranging from a 
2x2 grid, corresponding to 0.025 < q?, qL < 0.045 GeV/c , to 3 x 3, 4 x 4, etc. as shown in 
Figure 4. For each nxn grid, N = n 2 is the number of degrees of freedom and the standard 



S 



deviation is expected to be a = y/2/n. The strong dependence of the number of standard 
deviations from unity as a function of the range of the analysis is brought out clearly in 
Figure 4. 




FIG. 4. Number of standard deviations of (x 2 ) from unity for increasing number of 
bins. In part (a), and At were optimized, whereas in (b), only Rt was. 

Although less striking than in the pion case, we see from Figure 4 that it still is possible 
to separate the two scenarios, although none of them could be considered as a very good fit. 
This situation could become better in the near future with improving statistics, which could 
allow for smaller bin sizes. It is clear, however, that the non-resonance picture is closer to 
the preliminary data in all the range significant for interferometry, i.e., in the smaller domain 
Qt,Ql < 0.20 GeV/c, where the correlation function deviates significantly from unity. The 
two models yield similar fits in terms of x 2 f° r Imax > 200 MeV/c because in that large 
domain both models trivially predict nearly unit correlation functions. 

In Figure 4.(b) two curves signaled with Ar = can be seen. This result corresponds 
to fixing the decoupling time interval to zero (instantaneous freeze-out) and searching for 
the optimized value of Rt- This test was performed following the suggestion in Ref. [§, 
according to which the preliminary 3-D experimental analysis in q , q?, q^ returned results 
for Rt, Rt compatible with values obtained by the 2-D analysis, although the value of Ar 
found could either be Ar ~ Rt or Ar = 0! We should recall that the experimental analysis 
had to project data in large bins (width of 180 MeV) 0, in order to have enough statistics 
in the time direction which, by itself, would be responsible for dramatically weakening the 
interferometric signal. However, that ambiguous conclusion regarding the time interval was 
reached when a static Gaussian space-time parameterization was used in the experimental 
fit. We then decided to test what would be the response of the method to it, since we consider 
a different class of models, in which the longitudinal expansion is taken into account. The 
region of Rt searched in its optimization, keeping Ar = 0, is shown in Figures 2 (a) and 
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(c). The optimized values for Rt are shown in Table 1. The errors, estimated with the aid 
of the asymptotic form of the P{x 2 ) distribution, would not apply to this case. The reason 
is that, for considering the average x 2 over N bins as a random variable with unit mean and 
rms width o = y2/iV, as discussed above, the x(hj) should be normal random variables 
with zero mean and unit rms width. We estimated this distribution in each case by running 
CERES for the optimized values. When fixing At to be zero, however, the assumption made 
about the x(hj) distribution was not verified. For this reason, we prefer to simply show the 
optimized values of Rt as approximate ones. The corresponding results can be seen both 
in Figures 3.(d) and (e) (with corresponding x 2 m (h) an d (i)), as well as in Figure 4.b. 
From this last one, we see that our model completely excludes the instataneous emission. 
In particular, even in the region where no correlation is expected (roughly for qx, qh > 200 
MeV), the deviation with respect to data continues to be enormous. 

Furthermore, the above analysis nicely illustrates the important and well-known @-|8| 
influence of the time spread in the effective transverse radius, Rt- Although its influence 
would be noticeable even for a static Gaussian parameterization of the space-time decoupling 
geometry ||, models considering expanding systems strengthens the effect HQ]. In the 
present analysis, the time influence on Rt can be inferred from the fact that the optimized 
radius increases considerably, trying to compensate for the strong constraint of zero emission 
time interval. For instance, when including the K* contribution, its finite lifetime tries 
to circumvent the problem by introducing a non-zero time spread through the resonance 
decay, albeit the optimized Rt is about twice the value without that constraint. This effect 
is, however, more dramatic in the non-resonance case, where no clear evidence about the 
location of the optimized value of Rt can be seen from Figure 2. (a), since there is no way 
out to accommodate the instant emission constraint. 

We conclude that the multi-dimensional analysis proposed in Ref. [0J, which has high 
resolving power in the domain of physical interest in the case of pion interferometry, can still 
be applied to the case of kaon interferometry, although with less resolving power, due to the 
absence of contribution from long lived resonances. Finally, the above two-dimensional x 2 
analysis indicates that, as far as the preliminary E859 data is concerned, expanding sources 
should be considered at the AGS energy range, since expansion enhances the influence of 
the emission time interval on the transverse dimensions of the source and, from the present 
analysis, kaon sources emitting instantaneously are discarded. This should also be considered 
as an alert against the common practice of employing the static Gaussian parameterization 
to fit data since, by using it, the interpretation of the corresponding extracted parameters 
could be misleading or even wrong. 
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